September 21, 2007 


Infrared and X-Ray Evidence for Circumstellar Grain Destruction 
by the Blast Wave of Supernova 1987A 

Eli Dwek 1 , Richard G. Arendt 2 , Patrice Bouchet 3 , David N. Burrows 4 , Peter Challis 5 , I. 

John Danziger 6 , James M. De Buizer 7 , Robert D. Gehrz®, Robert P. Kirshner 5 , Richard 
McCray 9 , Sangwook Park 4 , Elisha Polomski 8 , and Charles Woodward 8 

ABSTRACT 

Multiwavelength observations of supernova remnant (SNR) 1987 A show that 
its morphology and luminosity are rapidly changing at X-ray, optical, infrared, 
and radio wavelengths as the blast wave from the explosion expands into the 
circumstellar equatorial ring, produced by mass loss from the progenitor star. 

The observed infrared (IR) radiation arises from the interaction of dust grains 
that formed in mass outflow' with the soft X-ray emitting plasma component of 
the shocked gas. Spitzer IRS spectra at 5 - 30 pm taken on day 6190 since the 
explosion show that the emission arises from ~ 1.1 x 10 -6 M© of silicate grains 
radiating at a temperature of ~ 180± 29 K. Subsequent observations on day 7137 
show that the IR flux had increased by a factor of 2 while maintaining an almost 
identical spectral shape. Tire observed IR-to- X-ray flux ratio ( IRX ) is consistent 


1 Observational Cosmology Lab., Code 665; NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
U.S.A., e-mail: eli.dwekiiJnasa.gov 

2 Science Systems & Applications. Inc. (SSAI). Code 665. NASA Goddard Space Flight Center. Greenbelt 
MD. 20771. U.S.A. 

' DSM/ D APNI A/ Service d'Astrophysique. CEA/Saclay, F-91191 Gif-snr- Yvette; Patrice. Bouchet4Jcea.fr 

4 Department of Astronomy and Astrophysics, Pennsylvania State University. 525 Davev Laboratory. 
University Park. PA 16802. U.S.A. 

5 Harvard-Smithsonian, CfA. 60 Garden St... MS-19. Cambridge. MA 02138. U.S.A. 

f ’Osservatorio Astronomieo di Trieste, Via Tiepolo. 11. Trieste. Italy 

" Gemini Observatory. Southern Operations Center, c/o AURA. Casilla 603. La Serena, Chile 

^Department of Astronomy. University of Minnesota. Minneapolis, MN 55455. U.S.A. 

9 Joint. Inst.. Lab. Astrophysics. University of Colorado, Boulder, CO 80309-0440. U.S.A. 



- 2 ~ 

with that of a dusty plasma with standard LMC dust abundances. This flux ratio 
Inis decreased by a factor of ~ 2 between days 6190 and 7137, providing the first, 
direct observation of the ongoing destruction of dust in an expanding SN blast 
wave on dynamic time scales. Detailed models consistent with the observed dust 
temperature, the ionization fluence of the soft X-ray emission component, and 
the evolution of IRX suggest that the radiating silicate grains are immersed 
in a 3.5 x 10 6 K plasma with a density of (0.3 — 1) x 10 4 cm -3 , and have a 
size distribution that is confined to a narrow range of radii between 0.02 and 
0.2 fxm. Smaller grains may have been evaporated by the initial UV flash from 
the supernova. 

Subject headings: ISM: supernova remnants - ISM: individual (SNR 1987A) - 
ISM: interstellar dust - Infrared: general - X-rays: general 


1. INTRODUCTION 

On February 23, 1987, Supernova 1987A (SN 1987A), the brightest supernova since 
Kepler’s SN in 1604, exploded in the Large Magellanic Cloud (LMC). About ten years 
thereafter, the energy output from the supernova became dominated by the interaction of 
its blast wave with the inner equatorial ring (ER), a dense ring of gas located at a distance 
of about 0.7 lyr from the center of the explosion, believed to be produced by mass loss from 
the progenitor star. The ER is being repeatedly observed at optical wavelengths with the 
Hubble Space Telescope ( HST) (?), at X-ray energies with the Chandra X-ray Observatory 
(??), at radio frequencies with the Australian Telescope Compact Array (ATCA) (?), and in 
the mid-IR with the Gemini South observatory (??) and the Spitzer observatory (??). The 
morphological changes in its appearance in these different wavelength regimes (presented in 
Figure 6 of McCray 2007) reveal the progressive interaction of the SN blast wave with the 
ER . The interaction regions appear as hot spots in the HST images, representing the shocked 
regions of finger-like protrusions that were generated by Rayleigh- Taylor instabilities in the 
interaction of the wind from the progenitor star with the ER (see Figure 2 in McCray 2007). 

The X-ray emission is thermal emission from the very hot plasma, and consists of two 
main characteristic components: a hard (kT & 2 keV) component representing a fast shock 
propagating into a low density medium, and a soft (kT ss 0.3 keV) component representing a 
decelerated shock propagating into the denser protrusions (??). The optical emission arises 
from the gas that is shocked by the blast wave transmitted through the dense protrusions 
in the ER. ( ?) , and the radio emission is synchrotron radiation from electrons accelerated by 
the reverse shock (?). The mid-IR emission spectrum is comprised of line and continuum 
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emission. The lines most, probably originate from the optically bright dense knots. The 
continuum that dominates the spectrum is thermal emission from dust that, was formed in 
the post main sequence wind of the progenitor star before it exploded. This dust could either 
be located in the shocked X-ray emitting plasma and heated by electronic collisions or in the 
optical knots and radiatively heated by the shocks giving rise to the optical emission (??). 
Dust has also formed in the ejecta of the supernova about. 530 days after the explosion, as 
evidenced by optical and IR observations of SN1987A (????), but its current contribution 
to the total mid-IR emission is negligible (?). 

Because the Gemini 11.7 /xm image correlates well with both the X-ray and optical 
emission, the possibility that the dust is radiatively heated in the knots was considered in 
detail by ?. Estimated dust temperatures of ~ 125 K fell short, of the observed value of 
~ 180 K but given the uncertainties in the model parameters, this discrepancy could not 
firmly rule out this possible scenario for the location and heating mechanism of the dust. 
However, the combined IR and X-ray observations can be more readily explained if the dust, 
resides in the shocked regions of the finger-like protrusions that give rise to the observed 
soft X-ray emission (?). We will adopt this scenario as our working hypothesis, and derive 
a self-consistent model for the composition, abundance, and size distribution of the dust to 
explain the evolutionary changes in the IR and X-ray fluxes resulting from the collisional 
heating and the destruction of the dust grains by the ambient plasma. 

We first review the basic physical principles that determine the temperature of colli- 
sionally heated dust, and describe how the IR emission can be used to probe the physical 
condition of the X-ray emitting plasma. We also discuss what information can be derived 
from the comparison of the IR and X-ray fluxes from the gas (§2.1). The physics of dust 
particles in a hot gas is discussed in more detail by Dwek (1987) and Dwek & Arendt (1992). 
In §3 we present a simple analytical model for the evolution of the grain size distribution 
and total dust mass in the gas that is swept up by an expanding SN blast wave. In §4 
we present Spitzcr low resolution 5 — 30 //in IRS spectra obtained on days 6190 and 7137 
after the explosion. The IR spectra are used to derive the temperature and composition 
of the shock heated dust. We use IR and X-ray observations of the SN to constrain the 
grain size distribution and the time at which the SN blast wave first crashed into its dusty 
surroundings. The results of our paper are summarized in §5. 


2. The Infrared Diagnostics of a Dusty X-ray Plasma 

The morphological similarity between the X-ray and mid-IR images of SN 1987A sug- 
gests that, the IR emission arises from dust that is collisionally heated by the X-ray emitting 
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gas. Simple arguments presented below show that, under certain conditions, the IR luminos- 
ity and spectrum of a dusty plasma can be used as a diagnostic for the physical conditions of 
the gas and the details of the gas-grain interactions. Details of the arguments can be found 
in ? and ?. 


2.1. The Dust Temperature as a Diagnostic of Electron Density 

The collisional heating rate, 7-f(erg s -1 ), of a dust grain of radius a embedded in a hot 
plasma is given by: 

n^na 2 n :j Vj £j (1) 

3 

where nj is the number density of the j - th plasma constituent, Vj , its thermal velocity, and 
£j its thermally-averaged energy deposition in the dust. In all the following, we will assume 
that the ion and electron temperature are equal. Then v e » v icni) and the dust heating rate 
is dominated by electronic collisions. 

Let Efc v be the thermally- aver aged energy deposited by electrons in the solid. If most 
electrons are stopped in the dust then E^p is, on average, equal to the thermal energy of 
the electrons, that is, E^ ev oc T e , where T c is the electron temperature. On the other hand, 
if most incident electrons go entirely through the grains, then E^ v is proportional to the 
electron stopping power in the solid, defined as p~ l (dE / dx) . At the energies of interest 
here the electronic stopping pow r er has an energy dependence of ( dE/dx ) ~ E ~~ 1 ! 2 (?), or 
dE ~ E ~ 1 / 2 dx so that Ed ep ~ T a. 

The functional dependence of the grain heating rate on gas density and temperature is 
then given by: 

H ~ a 2 n c TfJ 2 when electrons are stopped in the grain (2) 

~ a 3 n e when electrons gothrough the grain 

1 '2 

where we used the fact that v e T (] ‘ . 

The radiative cooling rate, C , of the dust grain with temperature by IR emission is 
given by: 

£ = ir a 2 crTj ( Q ) (3) 

~ 7 TO 3 (JT^ 0 

where a is the Stefan-Boltzmann constant, and (Q) oc ciT 1 ] is the Planck-averaged value of 
the dust einissivity, Q( A) a \~ 0 , where the value of emissivity index, (3, is « 1 — 2. 
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In equilibrium, H = C, and the dust temperature dependence on plasma density and 
temperature can be written as: 

T,i ~ (— ) Tp' 2 when electrons are stopped in the grain (4) 

~ nj when electrons gothrough the grain 

where 7 = 1/(4 4 - (3). These simple arguments show that when the gas temperature is 
sufficiently high, and the grain size sufficiently small so that most electrons go through the 
grain, the dust temperature depends only on the plasma density. 
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Fig. 1. — Contour plot of the equilibrium temperature of 0.01 pin (left) and 0.10 pm (right) silicate grains 
as a function of electron density and temperature. Above temperatures of ~ 5 x 10 c K the 0.01 pm grains 
become transparent to the incident electrons, and the dust temperature is only a function of electron density. 
The larger 0.1 pm grains become transparent only at electron temperatures above ~ 3 x 10 7 K. We note 
here that a similar figure (Figure 15 in Bouchot et al. 2006) was mislabeled, and actually corresponds to 
contour levels of silicate dust temperature for grain radius of a = 0.0030 pm. 

Figure ?? depicts contour levels of the dust temperature as a function of electron density 
and temperature for 0.01 gm and 0.10 gm silicate grains. The figure shows that above a 
certain gas temperature, its value depending on the grain radius, most of the electrons go 
through the grain and the dust temperature is essentially determined by the electron density. 
Under these conditions, the IR spectrum from the collisionally-heated grains becomes an 
excellent diagnostic of the density of the X-ray emitting gas. 

2.2. The Stochastic Heating of Grains by Electronic Collisions 

When dust grains are sufficiently small, a single electronic collision can deposit an 
amount of energy in the dust that is significantly larger that its enthalpy, causing a surge in 
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Fig. 2. — The stochastic heating of silicate grains in a hot X-ray emitting gas characterized by a temperature 
of T (J = 10 (i K. and electron density n r = 1 cm -3 for dust grains of different radii. Left column: The 
temperature fluctuations as a function of time. Right column: The histogram of the fluctuations. As 
the grain size increases, the fluctuations get smaller, and the probability distribution of dust temperatures 
becomes strongly peaked around the equilibrium temperature of ~ 38 K. 
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dust temperature. If additionally, the time interval between successive electronic collisions is 
larger than the dust cooling time, the grain temperature will be fluctuating with time (??). 
Figure ?? depicts a simulation of the stochastic heating of 0.003, 0.01, and 0.03 fj,m silicate 
grains immersed in a hot X-ray emitting gas characterized a temperature T g = 10 6 K, and 
an electron density n e = 1 cm -3 . The left column shows the temperature fluctuations as 
a function of time, and the right column the histogram of the grain temperature. As the 
grain size increases, the fluctuations get smaller, and the histogram becomes strongly peaked 
around the equilibrium dust temperature of ~ 38 K, in this example. 
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Fig. 3. — Left panel: The cooling function of a dusty plasma via atomic processes (thick solid line) and 
via gas grain collisions. Calculations were performed assuming a single-sized population of 0.02. 0.1 //m and 
0.20 /mi grains (solid line) with a dust-to-gas mass ratio of 0.0062. which is the value in the local ISM of 
the bare silicate-graphite+PAH dust model of ?. The gas cooling rate per unit volume for both processes is 
given by L = ri* A (T). Right panel: The value of IRX for single-sized dust populations with radii of 0.02, 
0.1, and 0.2 //in with the same dust-to-gas mass ratio as the figure on the left. 


2.3. The Infrared to X-ray Flux Ratio 

Another important diagnostic of a dusty plasma is IRX, defined as the ratio of the IR 
to X-ray fluxes emitted by the gas (Dwek et al. 1987). If the dust is collisionally- heated 
by the gas then the total IR. flux, F/^, emitted from a gas volume V is proportional to 
n ( Ud A d{T g ) V , where rid is the number density of dust particles, and A d{T g ) is the cooling 
function (units of erg cm 3 s -1 ) of the gas via gas grain collisions. The total X-ray flux, Fy , 
from the same volume is proportional to n * A g (T (J ) V , where A 9 (T 5 ) is the cooling function 




of the gas via atomic processes. Thus 
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IRX = 



A d{T g ) 
A g (T 9 ) 


(5) 


Both cooling functions represent the energy losses through collisional processes, characterized 
by (a v E ) summed over all interactions in the plasma, where a is the cross section, v is the 
relative velocity of colliding species, and E is the energy lost in the process. 


For a given dust-to-gas mass ratio, that is, a fixed ( rid/n e ) ratio, IRX depends only 
on plasma temperature. Figure ?? (left panel) shows the behavior of the atomic cooling 
function of a gas of solar composition as a function of gas temperature. Also shown in the 
figure is the gas cooling function via gas-grain collisions for a gas with a dust-to-gas atom 
mass ratio Zj = 0.0062, and single-sized dust populations with radii of 0.02, 0.1, and 0.2 pm. 
The right panel of the figure presents the value of IRX for the same conditions. The figure 
shows that for soft X-rays {kT„ ~ 0.3 keV, T e ~ 3.5 x 10 6 K) this ratio varies between ~ 
3 and 20, depending on grain size. Each plasma temperature will have a different range of 
values, depending on the grain size distribution. Any deviation from these values will suggest 
that Zd is either depleted or overabundant with respect to the reference value adopted in 
the calculations. 


3. The Evolution of the Grain Size Distribution and Dust Mass 

Consider the propagation of a shock into a dusty medium with a constant number 
density, n 0 , and a fixed dust-to-gas mass ratio Z%. The shocked gas can be regarded as a 
reservoir that is continuously being filled with gas and pristine dust by the expanding blast 
wave. If dust grains were not destroyed, the postshock gas would maintain a constant value 
as the mass of shocked dust and gas evolve proportionally in time, with a functional 
dependence that depends on the geometry of the medium into which the blast wave is 
expanding. In the cast; where dust particles are destroyed by sputtering, the grain size 
distribution, and the dust-to-gas mass ratio in the shocked gas will evolve with time. 


3.1. General Equations 


A grain of radius ao that is swept up by the shock at some time to will at time t be 
eroded to a radius a(f) given by 

“ 55 a(<) = a « + i (!) M 


(6) 



where 


da 

dt 


r 


( 7 ) 
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y ; n j v j < ~ o 


where Yj is the thermally-averaged sputtering yield of the dust by the j-th gas constituent. 



Fig. 4. — The absolute value of the sputtering rate of silicate (Mg 2 Si0 4 ) dust grains moving through a hot 
gas of solar composition as a function of gas temperature. The curves are marked by the grain velocity (in 
km s _1 ). 


Dust particles swept up by a high velocity shock will move ballistically through the shock 
front and acquire a velocity relative to the shocked gas. The sputtering yield needs then 
to be averaged over a Maxwellian distribution of velocities that is displaced by the relative 
gas-grain motion from its origin in velocity space (?). Figure ?? shows the temperature 
dependence of the sputtering rate of silicate dust grains, calculated using sputtering yield 
parameters given by Nozawa et. al. (2006), moving with velocity v gr — 0, 50, 100, 200, and 
500 km s -1 through a hot gas of solar composition. In contrast to the heating of grains, 
their erosion by thermal and kinetic sputtering is entirely done by the ionic constituents of 
the gas. For dust grains with velocities > 500 km s" 1 and gas temperatures above ~ 10 6 K 
the sputtering rate is approximately constant and given by: 



—5 x 10 6 n//(cm 3 ) 

— 1.4 x 10 -4 n H (cnr 3 ) 


fxm yr 

A d - 1 


( 8 ) 


If the shocked gas maintains a constant composition and density then a dust grain of initial 
radius a 0 that is swept up by the shock at some time t' wall at time t have a radius a given 
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by: 

f' da\ , .. 

a — cio + I — I [t — t ) 
\at J 


(9) 


Equation (??) can be written in dimensionless form as: 

f = fo ~ = €o-(r}~ yf) 

•max 

where £ = a/a max , Co = a 0 /a max , V = t/r max , rjf = t f /r max , and 

^ max = ®mai | da/ dt | 


(10) 
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is the sputtering lifetime of the largest grain in the injected size distribution, which (using 
eq. ??) is numerically given by: 


■Tmax 


2 x 10 5 
7 x 10 3 


®maa:(M^) 
n# (cm -3 ) 

Gmax (A) 
n[{(cm ~' c *) 


yr 

d 


( 12 ) 


Let rid be the total number density of dust grains in the preshocked gas, and n^(a 0 ) da 0l 
the number density of grains with radii between a 0 and a 0 + da Q . We assume that the grains 
have a size distribution in the pre-shocked gas given by: 

J r 00 roo 

' «d(«o) dao = rid / f{a 0 )da Q (13) 

o Jo 

where f(a o) is the normalized size distribution. If the grain size distribution extends over a 
limited range of radii, a min < a 0 < a max , then f(a 0 ) = 0 for any a 0 < a min or a 0 > a max . 

Dust grains are continuously injected into the shocked gas by the expanding SN blast 
wave. The total number of shocked grains with radii a in the {a, a + da} interval at time t, 
Nd{a,t)da , is equal to the number of all dust particles of initial radius ao that were swept 
up at time t' (0 < < t) and sputtered during the time interval t' — t to radius a given by 

eq. (??). If V(t) is the growth rate of the volume of the shocked gas, then iV (/ (a. t) can be 
written as: f 

N d {a,t) = n d [ V(t')f(a 0 )dt' (14) 

Jo 

The lower limit of the integral, t = 0. corresponds to the time when the blast wave first 
encounters the dusty medium. 



The total mass of shocked dust at any given time t is given by: 


M d {t) = / 

** tt mt n 


\dafdt\t 


m d (a ) N d (a, t) da 


where md{a) = 47rpa 3 /3 is the mass of a dust grain of radius a. 


Equation (??) can be written in dimensionless form: 


NS, r\) = n d f V (rf) /(£ + r? - rf) dr]' 

Jo 


This integral is a convolution of the form: V(rf) * /(rjo — rj'), which can be numerically 
evaluated for arbitrary functions using Fourier transforms. 


3.2. A Simple Analytical Solution 


An analytical solution can be derived for a pre-shocked grain size distribution given by 
a power law in grain radius, and a power law time dependence of V. We write the grain size 
distribution as: 


/(a 0 ) = Ca 0 ‘ 


^min — a 0 ^ (l,nax 

otherwise 


where C = (k — l)/(a„,in 1 ~ “max 1 ) i s the normalization constant. 
The time dependence of V(t') can be written as: 


v{t’) = u 0 


where Vq is a proportionality constant, and a = 2 for a spherical blast wave expanding 
into a uniform interstellar medium (ISM), and a = 0 if the blast wave expands into a 
one-dimensional “finger-like” protrusion. 

The total number density of grains in the {a,a + da} radius interval is then given by: 

, r ( f v 


N d {a,t) = n d V 0 f (^—) f{a 0 )dt' 
Jo \ 7~inax J 


Using eq. (??) to change variables from t' to no, eq. (??) can be rewritten as: 


-I pa up 


N d (a,t ) - N d 


t a 

1 

^~max tirnax 


a 0 k dao 


(20) 
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where N d = rid Vo- 


The time dependence of N d (a, t) is contained in the limits on the integral over the grain 
size distribution. If the radius a is within the range of the injected grain size distribution, 
that is, Omin <a< a max , then a low = a, since only grains with radii larger than a could have 
contributed to N d (a,t). If the radius a is smaller than a nitn , then N d (a. t) is non-zero only 
if a + \da/dt\t exceeds a min , and a Uyw = a min . In other words, the most recent injection of 
grains that could have contributed to N d {a, t ) occurred at time t - At, where At is the time 
required to reduce the grain radius from a min to a. The largest grains that could have been 
sputtered to radius a during the time t is equal to a + \da/dt\t. However, the largest grain 
size cannot exceed a max • so the upper limit on the integral, a up , is determined by the smaller 
of these two quantities. To summarize: 


a low — max {a m m, o} 
dup = mill “f &max > d + 


dd 

dt 


( 21 ) 


For a shock expanding into a one-dimensional protrusion (a = 0) the solution to eq. 
(??) is given by: 


N d (a, f) a =o = N d 


da 

dt 


-l 


\n~ k+1 - a -fc+1 l 

l a law a up J 


( 22 ) 


(k- 1 ) 

At early times, when t < a\da/dt\~ l <C T max , a up — a , and ~ a, and the solution to eq. 
(??) becomes: 

N d (a,t) a ^ = N d tC a~ k . (23) 

At late times, when t > r max , a up = a max , and the asymptotic solution of (??) is: 


Nd ( d , t Tjjiqx )q — 0 A d 

= Nd 


da 

dt 

da 

dt 


(k-l) 




~k -\- 1 fc+11 


a > a n 


(24) 


-l 


C 


(k - 1) 


'ajn 1 - a ,nat l ] = constant a < a min 


For a shock wave expanding into a homogeneous medium (a= 2) the solution is given 


by: 


N d {a, t ) 


q-2 


N d 


da 

dt 


-l 


C < - 


+ 


[k 1) \ T inax Cl 


\ 2 

— ) (« 

max / 


A:+l — At+1 


low 


V‘) ( 25 ) 


t a 

+ 


+ 


{k 2) \T ma j) ci 

i 

(k - 3) 


-A: +2 


x low 


Up 


an 


-A.M-3 ___ n -k+3 
a low U up 
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Fig. 5. — Evolution of the grain size distribution with time, measured in units of r nmx . the sputtering 
lifetime of the largest grain with radius a m(tT in the size distribution. Calculations were performed for a 
grain size distribution characterized by an ~ a -3,5 power law in grain radii between 10 and 100 A. The grain 
destruction rate, da/dt , was taken to be 0.14 A d~h Bold lines are labeled by t/r nia x • Left: A spherical 
blast wave expanding into a one-dimensional protrusion (a — 0). Right: A spherical blast wave expanding 
into a uniform ISM (a = 2). 

At late times, when t » r max , the first term dominates, and the asymptotic solution of eq. 
(??) increases with time as t 2 : 

Nd(a,t > Trnaar)a= 2 ~ 

where 

0(a, t) = [a" fc+1 - a"^ 1 ] a > a min (27) 

^ max j COTlStCLTlt d ^ a>min 


da 

c 

^ i 

f * ) 

dt 

(k - 1) 1 

\ Tmax J 


G(a,t) (26) 


Figure ?? depicts the grain size distribution for different epochs. Select epochs, labeled 
by the dimensionless quantity t/T max , are represented by bold lines. Calculations were per- 
formed for an initial grain size distribution characterized by an ~ a -3 - 5 power law in grain 
radius between 10 and 100 A. The grain destruction rate, \da/dt\, was taken to be 0.14 A 
d“\ calculated for the sputtering rate of silicate grains in a hot gas with a temperature and 
density of ~ 10 6-8 K, and 1000 cm -3 , respectively 

The figure illustrates the dependence of the evolution of the gi'ain size distribution on 
the geometry of the ISM into which the blast wave is expanding. For a one-dimensional 
protrusion (a = 0), the figure (left) shows a clear convergence of the size distribution to a 
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fixed functional form and total number of grains for t/r max > 1. As the shock wave expands, 
the thickness of the shell of swept up dust increases with time. However, because of the 
finite grain lifetime, its thickness cannot exceed a value of AR s h ~ v s hT max , where v 8 h is the 
shock velocity. So the grain size distribution and total mass reaches a steady state limit. 
When the blast wave expands into a uniform medium, the shell of shocked dust reaches the 
same steady state thickness AR s h. However, since the surface of the shell increases as R 2 h , 
where R sfl is the radius of the blast wave, the mass of shocked gas will continue to increase. 
This is clearly depicted in the right panel of the figure, which shows that Nd(a,t) reaches a 
steady state, but continues to increase with time as t 2 . 

If grains were not sputtered in the shocked gas, then N®{a, i) da, the total number of 
dust grains in the {a, a + da} radius interval that are swept up by the shock at time t would 
be: 

C / t \ a+1 

N° d (a,t)da = N d r max — - — ( ) a~ k da (28) 

(a + 1) \T max J 

Their mass, is given by eq. (??) with N<j(a,t) replaced by the expression above, and 
with \da/dt\ set to zero. 



Fig. 6. — Evolution of the dust mass, M<*. the dust mass if grains were not sputtered, A/ , . and the fraction 
of the surviving dust. with time, measured in units of r max . the sputtering lifetime of the largest 

grain in the size distribution [r m « x (d) = 7.3 x 10 ' (i max (/mi) (cm -3 )] . Dust and gas parameters are 

identical to those used in Figure ??. The figure shows that the fractional change in M 4 /M® between two 
epochs constrains the grain size distribution and the density of the X-ray emitting plasma. Left: A spherical 
blast wave expanding into a one-dimensional protrusion (a = 0). Right: A spherical blast wave expanding 
into a uniform ISM (a = 2). 

Figure ?? shows the evolution of dust mass for a spherical blast wave expanding into a 
one-dimensional protrusion (q = 0; left), and into a uniform ISM (a = 2; right) as a function 
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of t/r nMx . As explained above, when a = 0, the mass of shocked dust, M d , reaches a constant 
limit for t » T max , whereas for a = 2 the mass of the shocked dust will increase as f 2 . If 
grains were not destroyed, the mass of swept up dust, M d would increase as t for a = 0, 
and as t 3 for a = 2. The figure also shows the evolution of the mass fraction of surviving 
dust grains, M d /M d . This mass fraction is proportional to the dust.-to-gas mass ratio in the 
shocked gas, and for a constant gas temperature and density, it is also proportional to IRX, 
the IR-to-X-ray flux ratio in the shocked gas. The figure shows that the fractional change 
in M d /M d between two epochs constrains the value of t 1imx given in eq. (??) which in turn 
depends on the grain size distribution and the density of the X-ray emitting plasma. As 
a reminder, T max {d) = 7000 a max ( A) /nH(cm~ 3 ). For example, given a plasma density, the 
value of T max will depend only on a max , the maximum grain radius. A small value of a max 
will imply a low value for T max , so that large changes in M d /M d occur over very short time 
scales. Conversely, large values of a mnx and r max , will cause changes in M d /M d to occur over 
very long time scales. 


4. Spitzer Infrared Observations of SNR 1987A 
4.1. The Evolution of the IR spectrum 

Figure ?? shows the 5 — 30 gm low resolution spectra of SNR 1987A taken on February 
4, 2004 (day 6190 since the explosion), and on September 8, 2006 (day 7137 since the 
explosion) with the Infrared Spectrograph (IRS) (??) on board the Spitzer Space Telescope 
(??). Analysis of the spectrum taken on day 6190 revealed that the IR emission originated 
from ~ 1.1 x 10 -6 Mq of silicate grains radiating at a temperature of 1801^5 K (Bouchet et 
al. 2006). These circumstellar grains were formed in the quiescent outflow of the progenitor 
star before it exploded. The total IR flux on day 6190 was 5.1 x 10 -12 erg cm -2 s" 1 (Bouchet 
et al. 2006), and increased after 947 days (day 7137)to 10.0 x 10" 12 erg crn~ 2 s" 1 . The right 
panel of figure ?? presents a comparison between the two spectra, both normalized to the 
same 10 fim intensity. The figure shows that the spectra are essentially identical, implying 
that the dust composition and temperature remained unchanged during the two observing 
periods. The lower curve in the figure is the ratio between the two spectra, emphasizing 
their similarity. The IR intensity increased by a factor of 2 between the two epochs. 
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Fig. 7. — Left: The evolution of the IR spectrum of SN1987A from April 2, 2004 (day 6190 since the 
explosion) to September 8, 2006 (day 7137) taken with the Spit.zer IRS (Bouchet et al. 2007; Arendt et 
al. 2007). Right: The smoothed spectra for days 6190 (solid line) and day 7137 (dashed line), normalized 
to the same brightness. The lower curve shows the ratio between the two spectra, with the horizontal line 
being the mean value. The figure shows that the dust spectrum increased by a factor of two between the two 
epoch, retaining essentially an identical spectrum corresponding to silicate grains radiating at an equilibrium 
temperature of 180t^ K. 


4.2. The IR-to- X-Ray Flux Ratio: Constraining the Dust Abundance in the 

ER 

A comparison of the IR and X-ray fluxes provides strong constraints on the dust abun- 
dance in the shocked gas. X-ray fluxes taken between days 6157 and 7271 with the Chan- 
dra, X-ray telescope (?) were interpolated for days 6190 and 7137 of the Spitzer observa- 
tions. The total X-ray flux on day 6190, corrected for an extinction H-column density of 
Nn = 2.35 x 10 21 cm -2 , is 2.1 x 10 -12 erg cm -2 s -1 , half of it radiated by the slow shock 
component (?). The IR emission originates from the slow shock component which is pene- 
trating the denser regions of the ER. This component comprises half of the observed X-ray 
flux. The resulting value of IRX on day 6190 is therefore 4.9 ± 1.1. 

The theoretical value for IRX in a gas with LMC ISM abundances, taken here to 
be 0.6 times solar (?), ranges from about 2 to 12 for soft X-rays with T e ~ 3 x 10 6 K. 
The dust abundance in the ER is therefore consistent with LMC abundances. Since the 
silicon abundance in the ER should not have been altered by stellar nucleosynthesis, this 
agreement suggests efficient condensation of silicate grains in the presupernova outflow. The 
dust abundance on day 6190 is a lower limit on the original pre-SN value, since some of the 
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dust may have been evaporated by the initial UV flash from the SN. 


4.3. The IR-to- X-Ray Flux Ratio: The Effect of Prior Grain Destruction by 

the Initial UV Flash 

HST images of the ER show that it is located at a distance of ~ 0.7 lyr (6.6 x 10 17 cm) 
from the SN. At this distance small dust particles can be evaporated by the initial UV flash 
that emanated from the SN (?). Their calculations suggest that silicate dust particles with 
radius less than ~ 0.02 /mi will be evaporated by the flash. A population of silicate grains 
with a a" 3,5 power law distribution in grain radius extending from 10 A to 0.2 /mi will loose 
about 30% of its mass. So had the flash not occurred, the value of IRX would have been 
~ 7. 


4.4. The Evolution of IR-to- X-Ray Flux Ratio: Evidence for Ongoing Grain 

Destruction by the SN Blast Wave 

The IR flux increased by only a factor of ~ 2 from day 6190 to 7137. In comparison, 
the extinction-corrected 0.50-2.0 keV flux increased by a factor of ~ 3 during the same time 
period to a value of ~ 6.4 x 10~ 12 erg cm -2 s" 1 (?). The fractional contribution of the soft 
X-ray component increased from 0.5 to 0.6, with no significant change in gas temperature 
( kT ~ 0.3 keV). All the increase in the soft X-ray flux can therefore be ascribed to an 
increase in the volume of the dense ( n e ~ 10 4 cm -3 ) component of the ER that was shocked 
by the SN blast wave. The evolution in the X-ray and IR fluxes and the resulting value of 
IRX are summarized in Table 1. A similar evolutionary trend was reported by ?? although 
absolute values of IRX differ from those reported here because of differences in the X-ray 
energy bandpasses and dust models used in the calculations. 

If grains were not destroyed, we would expect the IR intensity to increase by a similar 
factor. The smaller increase in the IR. flux, e.g. the decline in IRX, is a strong indicator that 
we are for the first time witnessing the actual destruction of dust in a shock on a dynamical 
timescale! If the dust composition and size distribution is uniform throughout the region 
of the ER that Inis been swept up by the shock, then IRX is directly proportional to the 
dust-to-gas mass ratio, Zj, or to the ratio between the actual mass of dust in the 

shocked gas, and the dust mass if grains were not destroyed. The magnitude of the decrease 
in IRX between the two epochs can therefore provide strong constraints on the grain size 
distribution in the preshocked gas of the ER, and the density of the X-ray emitting gas. 
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which determines the rate of grain destruction. 


4.5. Determining the Grain Size Distribution and Plasma Density 

If the preshocked grain sizes are too large, then the fractional mass of the dust that could 
be destroyed during the time interval of 947 days will be too small to account for the observed 
decrease in the value of IRX. Conversely, if the grain sizes were too small, most of the dust 
mass would be destroyed, giving rise to a significantly larger than observed decrease in IRX 
between the two epochs. The right combination of grain sizes and gas densities is therefore 
required to produce the observed decrease in IRX and dust temperature. Furthermore, the 
small range of dust temperatures ( Td « 165 — 200 K) implied from the mid-IR spectrum of 
the ER also constrains the grains to a rather narrow range of sizes. 

For a given gas temperature, the gas density is constrained by the size and temperature 
of the collisionally heated dust grains. The temperature of the soft X-ray component that 
heats the dust is about 3 x 10 6 K. If the dust grains are transparent to the incident electrons, 
then the dust temperature uniquely determines the gas density (see eq. ??). Inspection 
of Figure ?? suggests that 0.01 pm dust grains are transparent to electrons with mean 
temperatures above ~ 3 x 10 6 K. The observed dust temperatures of 180 K then requires 
the electron density to be about ~ 10 3 cm -3 . 

However, we do not know a priori if the dust grains are transparent to the incident 
electrons. If the grains are sufficiently large so that most of the incident electrons are 
stopped in the grains, then eq. (??) shows that for a given gas temperature, larger gas 
densities are required to obtain the same dust temperature that transparent grains would 
have. For example, Figure ?? (right) shows that gas densities of ~ 10 4 cm -3 are required to 
collisionally heat 0.1 pm dust grains to 180 K with a 3.5 x 10 6 K plasma. 

The degeneracy between the different combinations of grain radii and plasma densities 
required to heat the dust to ~ 180 K can be lifted by considering the additional observational 
constraints imposed by the evolution of IRX and the plasma ionization fluence derived from 
modeling the soft X-ray emission from the ER. 

Figure ?? depicts the evolution in the I RX (t 2 ) /IRX (t\) ratio as a function of t\ — to, 
where to is the time, since the explosion, when the SN blast wave first crashed into the 
dense material of the ER, and where — 6190 d and £2 = 7137 d correspond to the two 
epochs of near-simultaneous Spitzer and Chandra observations of the ER. The observed 
IRX{t 2 )/IRX{t\) ratio is ~ 0.53 ±0.16 (see Table 1), and shown as a horizontal dashed line 
in the figure. The grain size distribution used in the calculations is characterized by a a ~ 3 5 
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power law in grain radii, extending from a minimum grain size of 10 A to a value of a max 
of 0.01, 0.04, and 0.1 /mi. Results are presented for a SN blast wave expanding into a one 
dimensional protrusion (a = 0, left column) and a uniform medium (a — 2, right column) 
with densities of 10 3 cm -3 (top row) and 10 4 cm -3 (bottom row). 

The figure shows that the I RX (t 2 ) / / RX (ti) ratio attains its lowest value when t\ — to 
is small, that is, when the first, encounter of the ER with the SN blast wave occured just 
before ti, the first epoch of Spitzer observations. Since t\ is very close to to, very little 
grain destruction could have taken place during the t\ — to epoch. The value of IRX(t{) is 
therefore close to its pre-shock value. Consequently, any subsequent destruction would lead 
to relatively great changes in IRX at t = t 2 - Conversely, t\ — to attains its largest value 
when to = 0 (a physical impossibility because of the finite time required for the SN blast 
wave to reach the ER). Then, the relative change in IRX between the two epochs will be 
the smallest, and IRX^)/ IRX{ti) — ♦ 1. 

Observationally, we can associate t 0 with the first encounter of the SN blast wave with 
the appearance of the first hot spot in the HST image from April 1997, about 3700 days after 
the explosion (?). The soft X-ray light curve shows that the rise could have occurred between 
days 3700 and 6000. The first epoch corresponds to the first appearance of the optical knots, 
and the latter epoch corresponds to the time when the flux from the soft X-ray component. 
(kT ~ 0.3 keV) exceeded that from the hard component (kT ~ 2 keV) (?). From the mid-IR 
light curves (?), the energy output from the SN became ER dominated around day 4000. 
Adopting days 4000 to 6000 as a reasonable estimate for to gives a range of possible values 
of ti - t 0 » 200 - 2200 d. 

Inspection of Figure ?? illustrates the difficulty in constraining the grain size distribution 
or the astrophysical scenario (1-D protrusion or uniform ISM) because of the large uncer- 
tainty in the IRX U 2 ) / 1 RX (ti) ratio, which ranges from 0.37 to 0.69, and the large range in 
t\ — to- If we adopt the nominal value of I RX{t 2 )/ IRX (t\) = 0.53, then the models clearly 
favor small grain sizes if n gas = 10 3 cur 3 , regardless of the value of a, but are degenerate 
in the grain size distribution when n gas — 10 4 cm -3 . The evolution the I RX (t 2 )/ IRX (7 1 ) 
ratio alone is therefore not sufficient to lift the degeneracy in the different combinations of 
{n e , a} required to heat the dust to ~ 180 K. 

However, the value of to also determines the ionization age of the hot gas. The ionization 
fluence derived from modeling the soft X-ray spectra taken on days 6914, 7095, and 7271 
is given by n e t e > 10 7 cm -3 d (?). Taking to to be between days 4000 and 6000 gives an 
ionization time t e fa 1000 — 3000 d, and corresponding electron densities n e fa (10 — 3) x 
1() 3 cm -3 . These densities are high enough so that even small dust grains with radii ~ 10 A 
will be collisionally heated to their equilibrium dust temperature. Furthermore, at these 
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high densities, an equilibrium temperature of ~ 180 K can only be reached if the soft X-ray 
electrons are stopped in the grains (see Figure ??). The narrow range of grain temperatures 
then suggests that the grain size distribution should have a narrow range as well, since Td ~ 
a~ 7 (see eq. ??). The combined constraints on the gas density (n e ta (0.3 - 1) x 10 4 cm -3 ), 
gas temperature (T e « 3 x 10 6 K), and dust temperature (Td ~ 165 — 200 K) can therefore 
be used to determine the range of grain sizes in the ER. From contour plots similar to Fig. 1, 
we find that the grain radii can be as small as 0.02 /am if the gas density is 3 x 10 3 cm -3 and 
as large as 0.2 /urn if the gas density is 1 x 10 4 cm -3 . Dust temperatures for these two cases 
are 200 and 165 K, respectively. 

A model consistent with the observed dust temperature, the ionization fluence of the 
soft X-ray emission component, and the evolution of the IRX ratio is depicted by the 
red curves in figure ??. They depict the evolution in the IRXfa)/ IRX(t\) ratio for a 
population of silicate grains with a narrow a~ 3 5 power law distribution in grain radius with 
{a m)n , a max } = {0.08, 0.1} /um, in a gas with a density of 10 4 cm -3 . Taking the uncertainty 
in the IRXfe)/ IRX(t{) ratio into account gives a value of t 0 « 4200 — 5800 d for a = 0, 
consistent with X-ray and mid-IR observations. An lower limit of t 0 4200 d is inferred for 
the same dust model when the blast wave is propagating into a homogeneous medium. 


4.6. Implications for Determining the Mass of the Circumstellar 
Environment of SN1987A Using Light Echoes 

It is interesting to compare the dust properties derived for the ER with those derived 
for the progenitor’s circumstellar environment from studies of the evolution and intensity of 
light echoes created by the scattering of the optical light from the supernova by the dust 
grains. At any given time, all points with equal delay time lie on an ellipsoid of revolution 
with the SN at one focal point and the observer at the other. Unfortunately, the ellipsoid at 
the earliest epoch at which the echoes were observed was outside the ER (set' Figure lln in 
?). As a result, the light echoes probed only the circumstellar and interstellar media exterior 
to the ER. 

Assuming cylindrical and reflection symmetry, ?? derived a model for the morphology of 
the scattering medium consisting of: (1) a peanut-shaped contact discontinuity (CD) between 
the red supergiant and main- sequence winds from the progenitor star; (2) a structure called 
Napoleon’s Hat (NH) constituting the waist of this peanut; and (3) the two outer rings of 
the circumstellar shell (CS) that define the hourglass that is pinched by the ER. To model 
the scattered light ?? used the ? model for interstellar LMC dust with grain radii ranging 
from an upper limit of 0. 2-2.0 /am to a lower limit of 0.00035 /am. By varying the relative 
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silieate-to-carbon dust, mass ratio while maintaining an LMC dust-to-gas mass ratio that is 

0.3 times the value of the local interstellar medium, they estimated a total nebular mass of 
1.7 M & . They also found that the gas density increases, the maximum grain size decreases, 
and the silicate-to-carbon dust mass ratio increases as the echo samples material that is 
closer to the SN. The ER, with its population of smaller pure silicate grains, is consistent 
with this trend. The higher value of the minimum grain size in the ER may be the result of 
its proximity to the SN which caused the evaporation of grains smaller than 0.02 pm by the 
initial UV flash. Finally, the dust abundance in the ER is consistent with that adopted by 
? for the nebula, supporting their derived value for the nebular mass. 

5. Summary 

The interaction of the SN 1987 A blast wave with the complex structure of the ER. has 
given rise to rapid evolutionary changes in the X-ray, optical and mid-IR morphology of 
the emission. The Gemini-S mid-IR images have established that the IR emission originates 
from dust in the ER. that is swept up by the SN blast wave, and collisionaly heated by a 
soft X-ray component which has a temperature of 3.5 10 6 K, and an ionization fluence of 
n e t > 10' cm -3 d. The Spitzer infrared observations provide important complementary 
information on the evolution of the interaction of the SN blast wave with the ER and the 
properties of the dust in the hot X-ray emitting gas. The results of our analysis are as 
follows: 

1. Spitzer spectral observations on day 6190 after the explosion revealed that the dust 
consists of silicate dust grains radiating at an equilibrium temperature of ~ lSOif® K. 
Subsequent observations on day 7137 revealed that the IR flux increased by a factor 
of ~ 2, with the same dust composition and temperature remaining the same (Figure 
??). 

2. The dust grains could attain this temperature if they are sufficiently small (a ~ 
0.01 pm) to be transparent to the incident electrons and are immersed in a relatively 
low density plasma with n e ~ 10 3 cm -3 . Alternatively, the dust grains could be large 
enough to stop the incident electrons (a ~ 0.1 pm) if they are immersed in a higher 
density plasma with n e ~ 10 4 cm -3 . 

3. The degeneracy in the {n e , a} combinations required to heat the dust to ~ 180 I< can 
be lifted by considering the additional constraints imposed by the observed evolution 

of the infrared- to- X-ray flux ratio ( IRX ) and by the ionization fluence obtained from 
modeling the soft X-ray emission component . 
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4. The value of IRX decreased by a factor of ~ 0.54 between days 6190 and 7137, sug- 
gesting that we are witnessing the effects of grain destruction on a dynamical timescale 
of the remnant. This decrease can be used to constrain the grain size distribution, time 
of impact of the SN blast wave with the ER, and the density of the X-ray emitting gas. 

5. To follow the evolution of IRX, we constructed a model for the evolution of the grain 
size distribution in the shocked gas. In the model, pristine dust is continuously injected 
into the hot gas by the expanding SN blast wave, and destroyed by thermal and kinetic 
sputtering behind the shock. The evolution of the grain size distribution resulting 
from the combined effect of dust injection and destruction is presented in Figure ?? 
for different geometries of the medium into which the blast wave is expanding. 

6. The evolution in IRX represents the changes in the dust-to-gas mass ratio in the 
shocked gas resulting from grain destruction (see Figures ?? and ??). A given decrease 
in IRX between two epochs constrains the epoch at which the dust is first swept up 
by the SN blast wave (which determines the ionization time), the plasma density, and 
the grain size distribution. Given a plasma density, if the grains are too small they 
will be rapidly destroyed, giving rise to large changes in IRX on short time scales. 
Conversely, if the grains are too large, the interval of time over which IRX drops by 
the observed factor will become too long. 

7. A self-consistent picture that emerges from the application of the model to the com- 
bined X-ray and IR observations is that of a SN blast wave expanding into a dusty 
finger-like protrusion of the ER. The dust in the preshocked gas consists of pure silicate 
dust with a normal LMC dust-to-gas mass ratio and a grain size distribution limited to 
radii between 0.02 and 0.2 /ttm, sufficiently large to stop the incident electrons. Smaller 
grain sizes may have formed in the mass outflow from the progenitor star but were 
probably vaporized by the initial UV flash from the SN. The SN blast wave crashed 
into the ER between days 3600 and 6000 after the explosion giving rise to the observed 
soft X-ray emission. Typical temperatures and densities of the soft X-ray emitting gas 
are ~ 3 x 10 6 K and (0.3 — 1) x 10 4 cm -3 , consistent with those required to collisionally 
heat the dust to a temperature of ~ 180 K. The plasma parameters and grain size 
distribution are consistent with the amount of grain destruction needed to account for 
the observed decrease in the IRX flux ratio between days 6190 and 7137. At these 
gas densities, the onset of grain destruction occurred about 1200 - 2000 days before 
the first Spitzer observations, consistent with the rise in the soft X-ray flux and the 
ionization time derived from X-ray models. The absolute value of IRX is consistent 
with that expected from a dusty plasma with a typical LMC dust-to-gas mass ratio, 
taken here to be 60% of the value of the local ISM. 
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Further Gemini and Spitzes observations of SNR 1987A are in progress which, with 
combined X-ray observations, will shed further light on the nature of the morphology and 
dust properties of the circumstellar medium around the SN. 
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Fig. 8. — Evolution of the ratio I RX fa) / 1 RX (t i) as a function of the time difference ti — to . The time 
to is the time since the explosion when the SN blast wave first encountered the dusty equatorial ring (ER). 
The times 1 1 and t .2 correspond to the two epochs of near-simultaneous Spitzer and Chandra observations 
of the ER. respectively, days 6190 and 7137 since the explosion. The dashed horizontal line depicts the 
nominal value of the I RX{ti) / 1 RX{t\) ratio which is 0.53 ± 0.16 (see Table 1). The curves are labeled by 
the maximum grain size (in A) of the distribution. The minimum grain size was taken to be 10 A in all 
cases. The red curves correspond to a narrow grain size distribution with {a mr „. a n ui.v} — {800. 1000} A. 
Figures are also labeled by the value of n gas . the density of the shocked gas. Left column: A spherical 
blast wave expanding into a one-dimensional protrusion (a — 0). Right column: A spherical blast wave 
expanding into a uniform ISM (a = 2). 




Table 1. Observed X-ray and Infrared Fluxes From SN 1987A 


day number 1 

X-ray flux 2 

fsoft 3 

IR flux 

IRX flux ratio 4 

6190 

(2.1 ±0.32) 10- 12 

0.50 

(5.1 ±0.9) 10- 12 

4.9 ±1.1 

7137 

(6.4 ±0.32) 10" 12 

0.60 

(1.0 ±0.18) 10~ n 

2.6 ±0.5 


1 Since the explosion 

2 X-ray flux in the 0.5-2. 0 keV band, interpolated to the epochs of the Spitzer 
observations and corrected for an extinction column density of Nh = 2.35 10 21 . 
The error represents the uncertainty in the pile up correction factor (Park et al. 
2007). 

3 The fraction of the 0. 5-2.0 keV flux that arises from the soft (kT ~ 0.3 keV) 
X-ray component (?). 

4 The ratio of the IR to soft X-ray flux from the SN. The IRX flux ratio has 
decreased by a factor of 0.53 ± 0.16 from day 6190 to day 7137. 



